Resonance-assisted tunneling in three degrees of freedom without discrete symmetry 
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We study dynamical tunneling in a near-integrable Hamiltonian with three degrees of freedom. 
The model Hamiltonian does not have any discrete symmetry. Despite this lack of symmetry we 
show that the mixing of near-degenerate quantum states is due to dynamical tunneling mediated 
by the nonlinear resonances in the classical phase space. Identifying the key resonances allows us 
to suppress the dynamical tunneling via the addition of weak counter-resonant terms. 
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Tunneling is a phenomenon that is forbidden by clas- 
sical mechanics but allowed by quantum mechanics. In 
general, any flow of quantum probability between (ap- 
proximately) equivalent yet classically disconnected re- 
gions constitutes tunneling. The classical regions could 
be disconnected due to barriers in coordinate space, mo- 
mentum space or, more generally, in the classical phase 
space. In the cases where tunneling occurs despite the 
absence of obvious energetic barriers it is called dynam- 
ical tunneling\^; the barriers now arise due to certain 
exact or approximate constants of the motion and hence 
are naturally identified in the underhringclassical phase 
space . Considerable theoretical |lL uL lift |£ la In LL la lit Li] 
and experimental^! works have established that tunnel- 
ing between quantum states localized on two symmetry- 
related regions in the phase space can be strongly in- 
fluenced by the classical stochasticity (chaos-assisted 
tunneling^) and/or by the intervening nonlinear res- 
onances (resonance- assisted tunneling 6] ) . In the for- 
mer case, phase space is mixed regular-chaotic and the 
splittings show marked dependence on the nature of the 
chaotic states which couple to the tunneling doublets 2, 
0,0]. In the latter case with near-integrable phase space, 
the splittings depend delicately on the various resonance 
islands bridging the degenerate states |a, |fj, la 111 Il0| . 
Clearly, a quantitative scmiclassical theory still elusive, 
requires one to identify key structures in the phase space 
on which the theory is to be based. In this regard there 
is increasing evidence 0, 0] that the classical nonlinear 
resonances might play a central role in near-integrable as 
well as mixed phase space situations. 

However, most of the studies thus far have been 
on two degrees of freedom (dof) systems with discrete 
symmetries Does the resonance- assisted tunneling 
viewpoint hold in systems with three or more dof which 
lack discrete symmetries? The main motivation for 
our study comes from suggestions 8] put forward in the 
molecular context - can dynamical tunneling provide 
a route for mixing between near-degenerate states and 
hence energy flow between regions supporting qualita- 



tively different types of motion? In addition, notwith- 
standing the difficulties associated with visualizing the 
multidimensional phase space, dynamics in three or more 
dof has features that cannot manifest in the systems stud- 
ied up until now|l3|. In this letter we attempt to un- 
derstand dynamical tunneling in a model nonsymmetric, 
near-integrable three dof system. We show that mixing 
of near-degenerate states occurs via dynamical tunneling 
mediated by nonlinear resonances and the mixing can be 
suppressed by adding weak counter-resonant terms. 
We study the Hamiltonian 



H = Hn 



- £ K mr [(4r M " (a 3 ) 7 " M 5 " + h.c], (1) 



describing four coupled modes j = 1, 2, 3, 4 with 



H = Y^^i n 3 + 



x jj n j , 



j<k 



(2) 



and -ffo| n ) = -E"nl n )- Although eq. ^has been inspired in 
the molecular context^], similar multiresonant Hamil- 
tonians arise in a variety of svstems|l5j. The occupa- 
tion number of the j th mode, rij = djdj, is expressed 



in terms of the harmonic creation 



and destruction 



(a.,) operators. The perturbations are characterized by 
m. r = (a r , — /3 r , — 7r, — o~r) with strengths K mr . The clas- 
sical limit of eq. generated via the correspondence 
Ij exp(i9j), is the following Hamiltonian: 



W(I,0) = W (I) 
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I?*ltl2 r Ii r cos(m r -0). 

(3) 

(I, 6) are the classical action-angle variables of 7io and 
hence the perturbations correspond to classical nonlin- 
ear resonances. The parameter e has been introduced 
for a perturbative analysis (see below). We restrict our- 
selves to three perturbations ni! = (1,-2,0,0), m 2 = 
(1,-1,-1,0), and m 3 = (1,-1,0,-1). This allows for 
a clear study of the role of the specific resonances in dy- 
namical tunneling. The existence of a conserved quan- 
tity P = ni + (ri2 + «3 + ri4)/2, with the classical analog 
P c = h + (h + h + h) /2, implies that the 4-mode system 
has effectively three dof. The eigenstates, eigenvalues, 
and the resulting mean level spacing of H are denoted 



2 



by | a), E a and AE respectively. In the units appropri- 
ate for the model Hamiltonian[l4J the Heisenberg time 
is given by t h = (2ncAE)^ 1 with c being the speed of 
light. 

We are interested in the fate of a set of near-degenerate 
zcroth-order states in the presence of weak perturbations, 
K m JAE = fc mr < 1. Consider states |n), |n'), |n"), . . . 
such that E^ w i=s E^,. «... with average energy 
E and E° € (E - AE/2,E + AE/2). Certain states, 
among the set of near-degenerate states, mix since they 
are directly connected to each other via one of the pertur- 
bations. The nonlinear resonances in eq.|3]do mediate the 
mixing via dynamical tunneling. However, in this work 
we will focus on states that are not directly coupled by 
the resonances in eq. [21 in order to show that even very 
weak, induced resonances can lead to substantial mixing 
that can be associated with dynamical tunneling. Quan- 
tum mechanically, the extent of mixing of a zeroth-order 
state |n) can be gauged by computing the survival prob- 
ability P n Q(t) and the inverse participation ratio (IPR) 

P nQ (t) = Knle-^ln}! 2 ^^^^^— (4) 

a, j3 

1 f T 

a n = lim - / P nQ {t)dt = YVL, (5) 

with p na = |(a|n)| 2 and uj a p = (E a - E@)/h. If a n < 1 
then |n) mixes extensively with other zeroth-order states. 

Specifically, we investigate a set of zeroth-order states 
around E ps 3542. 5A£ and P = 8. This choice of 
E is motivated by the existence of a number of near- 
degenerate states and qualitatively similar behavior is 
seen at different values of E as well. We select states 
that are not directly coupled by the perturbations in 
eq. ^ but nevertheless have IPR smaller than one. In 
Fig. ^ we show the survival probabilities for four such 
zeroth-order states |o) = |0,11,1,4), \b) = 0,11,2,3), 
|c) = |0,12,2,2), and \d) = |0,13,1,2) with cr n = 
0.97,0.51,0.40, and 0.74 respectively. The crucial thing 
to note is that |6),|c), and \d) mix amongst themselves 
over long times. The corresponding classical calculations, 
shown in Fig.^middle row, indicate long time trapping. 
Thus the observed mixing between the states is classically 
forbidden and corresponds to dynamical tunneling. In 
Fig.0the variation of the energy levels with the coupling 
parameter k mi is shown to indicate the lack of avoided 
crossings between the states of interest. Fig. also shows 
the spectral intensities p na and in every case we see two 
clumps of lines - one at the origin and another w 60AE 
away. A purely quantum explanation invokes the second 
clump of states, the virtual or off-resonance states, which 
provide a "vibrational superexchange" pathway for the 
mixing jlH|. The virtual states themselves have a n « 1 
and hence do not mix significantly. It is particularly 
striking to note that neither p nct nor the energy level vari- 
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FIG. 1: Quantum survival probability (top row) P n Q of the 
states |n) = |0, 11,1,4) = |o) (black), |0, 11,2,3) = \b) 
(red), |0,12,2,2) = |c) (green), and |0, 13, 1, 2) = |d) (blue). 
Time r is measured in units of the Heisenberg time and 
0.5K mi = K m2 = K m3 w 0.2AE. AE w 4.44 cm" 1 for 
eq.Qwith P = 8. The cross probabilities |(n'|n(t))| 2 are also 
shown. The middle row shows the classical analog, P n c, of 
PnQ- PnC indicates trapped motion and is qualitatively dif- 
ferent from P n Q. The bottom row shows the quantum P n g 
with the addition of weak counter-resonant terms (eq.|7J. The 
dynamical tunneling essentially shuts down, proving the im- 
portance of the induced resonances. 



ations suggest any differences between the zeroth-order 
states in contrast to the observations in Fig. ^ We now 
show that a relatively simple explanation can be given in 
terms of resonance-assisted tunneling on the energy shell. 

In the resonance-assisted tunneling scenario the mix- 
ing between, for example, \b) and |c) can be mediated 
by a 1:1 resonance involving modes j = 2 and 4 i.e., 
a resonance vector = (0,1,0,-1). The Hamilto- 
nian in eq. |3| does not have explicitly but it can be 
induced by mi and m 3 . Similarly, m^i = (0,0,1,-1) 
and m i2 = (0,1,-1,0) can be induced by the reso- 
nances in eq. [21 The resonances can be visualized by 
constructing the Arnol'd web^j| at E w E and fixed 
P, i.e., the intersection of the various resonance planes 
m r • dn (l)/dl = with the energy shell H (I) ~ E. 
For near-integrable systems the energy shell, resonance 
zones, and the location of the zeroth-order states can 
be projected onto a 2-dimensional space of two indepen- 
dent frequency ratios. The "static" Arnol'd web based on 
Ho highlights the various possible resonances and their 
topology on Hq(T) ~ E. However, from the tunneling 
perspective it is crucial to determine dynamically rele- 
vant part of the static web at E w E. This "dynamical 
web" is determined via a wavelet based local frequency 
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FIG. 2: Eigenvalues E = E a /AE versus varying coupling 
strength k mi = K m JAE. K m2 = K m3 « 0.2 AE. The 
circles correspond to the eigenstates having the largest con- 
tribution from the specific zeroth-order states (Fig. 0. The 
dotted lines show a similar calculation using eq. [7| indicat- 
ing no qualitative change in the level motions. The bot- 
tom panel shows the overlap intensities p na = |(a|n)| 2 versus 
Ae n = (E° — E a )/AE. Note the log-scale for the intensities 
and the cluster of lines in all the plots around Ae n = 0, and 
60. 



analysis (l7| of the Hamiltonian eq. [31 Briefly, initial con- 
ditions (I, 6) satisfying 9) ss E are generated and 
the trajectories are followed in the frequency ratio space 
(O1/O3, Oi/ O4) . The frequencies Qfc(t) are computed 
along the trajectory (1, 8)(t) by performing the wavelet 
transform of Zk(t) = y / 2Ik(t) ex-p(i8k(t)). The total num- 
ber of visits in a given region of the ratio space gives a 
density plot representing the web which highlights dy- 
namically significant regions at a given energy. The re- 
sulting "dynamical web" is shown in Fig. along with 
the location of the relevant zeroth-order states. Apart 
from highlighting the primary resonances mi, 1113 the fig- 
ure also indicates the existence of the induced resonances 
niii,mi2, and separating the states. The states are 
located close to the junction formed by the three induced 
resonances and far away from the primary resonances. 
Hence Fig. 01 supports the notion that \b), |c), and \d) 
are mixed due to dynamical tunneling mediated by the 
induced resonances. 

In order to conclusively establish the role of the in- 
duced resonances it is necessary to extract their strengths 
by perturbatively^ 0, Hof removing the primary reso- 
nances in eq.|2|to 0(e). As a result we obtain the effective 
Hamiltonian containing the induced resonances at 0(e 2 ) 
which are approximated by effective pendulums. For in- 
stance, one obtains the effective pendulum Hamiltonian 
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FIG. 3: The "dynamical" Arnol'd web at E = E generated by 
propagating 25000 classical trajectories for r « 40. The pri- 
mary resonances, mi (red), rri2 (black) and 1x13 (magenta), 
as predicted by Tto are superimposed for comparison, m.2 
and rri3 intersect on the energy shell giving rise to the in- 
duced resonance m^i = (0,0,1,-1). 1112 and 1H3 interact 
with mi to give rise to the two other induced resonances 
m 4 2 = (0,1,-1,0) and m« = (0,1,0,-1) respectively. The 
nearly degenerate zeroth-order states (filled circles) are lo- 
cated close to the induced resonances (arrows) which lead to 
state-mixing (cf. Fig. via resonance-assisted tunneling. 



-,(24) (A 24 - A 24 j nv „, 

He ff = 2M^ + V24 " 024 ' (6) 

appropriate for the induced resonance A mi3 with K24 ~ 
2/4 and 2^24 ~ (#2 — #4)- The resonance center is denoted 
by A^2 4 . The coupling V24 can be expressed in terms 
of the conserved quantities I^,P C , and P24 = I2 ~\~ Ia 
and the resulting tunneling time T tun « {2 j nAE/V24) 
agrees well with Fig. ^ The strengths estimated via a 
pendulum approximation can be translated back to ef- 
fective quantum strengths A m and our analysis reveals 
that the induced resonances are more than an order of 
magnitude smaller then that of the primary resonances 
(Amj, < |A mi2 | w |A mi3 | « 0.035A' mi ). Note that m l2 
and come with a negative sign as opposed to m^i . It 
is known; 6] that for significant mixing the states must lie 
symmetrically with respect to the center of the mediating 
resonance zone. Among the states considered, \b) and |c) 
satisfy the criterion very well and hence enhanced mixing 
between them is seen in Fig.^ The state \a) is not sym- 
metrically located with respect to \b) and thus, combined 
with the very small strength A mil , the induced resonance 
mji is ineffective. Now consider modifying eq. H accord- 
ing to 

H' = H+ |A mi2 1(4^3 + h.c.) + |A mi3 |(4a 4 + h.c), (7) 

where we have added terms to counter the induced res- 
onances. The reasoning is simple - if the induced reso- 
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FIG. 4: IPR (eq.0 of the near-degenerate states with respect 
to H (circles), H' with the correct sign (squares), and H' with 
the opposite sign (triangles). |Am i2 | ~ |A mi3 | ~ 0.014AZ? 
in eq. |7| Note the importance of the sign of the induced 
resonances and the shutdown of tunneling for the states of 
interest (filled symbols). Other nearby states (open symbols) 
are least affected by the counter-resonant terms. 



nances are truly mediating the dynamical tunneling then 
adding the counter-resonances should suppress the tun- 



neling. Moreover, since A 



A m ,, < K 



m i2 i ' ( m;3 



quantities 

like mean level spacing, eigenvalue variations (cf. Fig.EJ) 
and spectral intensities show little change as compared 
to the original system. Despite this, as shown in Fig. 
bottom row, the survival probabilities for \b), |c), and \d) 
indicate an almost complete shutdown of dynamical tun- 
neling. In Fig. 01 we show that the IPR also indicate 
the shutdown of tunneling whereas IPR of other nearby 
states do not change with the inclusion of the counter 
terms. The result of adding the counter terms with the 
same magnitude but opposite sign shows (cf. Fig. |3J en- 
hanced mixing and this establishes that the induced res- 
onances rrii2 and are responsible for the dynamical 
tunneling seen in Fig. ^ 

In conclusion, this work shows that significant mix- 
ing between near-degenerate states due to resonance- 
assisted tunneling can be expected in very general situa- 
tions. In addition, by suitable local modifications of the 
phase space, complete control of the dynamical tunneling 
can be attained. In light of a recent work the counter- 
resonances can be thought of as weak control terms 
and in nonautonomous systems this suggests the possi- 
bility of control via additional weak driving fields with 
particular attention to the relative phases between the 
fields 0| . The model system studied herein is certainly 
not in the deep semiclassical limit, perhaps reason enough 
to argue against competition from classical transport 
mechanisms, and yet the importance of the nonlinear res- 
onances is clear. Further work in the deep semiclassical 
regime, more closely approaching the molecular systems, 
is in progress. 
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